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ABSTRACT 

We present mult i- wavelength analysis around mid-infrared (MIR) bubble N14 to probe 
the signature of triggered star formation as well as the formation of new massive 
star(s) and/or cluster (s) on the borders of the bubble by the expansion of the Hll 
region. Spitzer-IRAC ratio maps reveal that the bubble is traced by the polycyclic 
aromatic hydrocarbon (PAH) emission following an almost circular morphology ex- 
cept in the south-west direction towards the low molecular density environment. The 
observational signatures of the collected molecular and cold dust material have been 
found around the bubble. We have detected 418 young stellar objects (YSOs) in the 
selected region around the bubble N14. Interestingly, the detected YSO clusters are 
associated with the collected molecular and cold dust material on the borders of the 
bubble. One of the clusters is found with deeply embedded intermediate mass and 
massive Class I YSOs associated with one of the dense dust clumps in the east of the 
bubble N14. We do not find a good agreement between the dynamical age of the Hll 
region and the fragmentation time of the accumulated molecular materials to explain 
possible "collect-and-collapse" process around the bubble N14. Therefore, we suggest 
the possibility of triggered star formation by compression of the pre-existing dense 
clumps by the shock wave and/or small scale Jeans gravitational instabilities in the 
collected materials. We have also investigated 5 young massive embedded protostars 
(8 to 10 Mq) and 15 intermediate mass (3 to 7 M ) Class I YSOs which are asso- 
ciated with the dust and molecular fragmented clumps at the borders of the bubble. 
We conclude that the expansion of the Hll region is also leading to the formation of 
these intermediate and massive Class I YSOs around the bubble N14. 

Key words: dust, extinction - Hll regions - ISM: bubbles - ISM: individual objects 
(IRAS 18134-1652) - stars: formation - stars: pre-main sequence 



1 INTRODUCTION 

Massive stars (M > 8 Mq) have the ability to interact 
with the surrounding cloud with their energetic wind, UV 



ionizing radiation and an expanding Hll region (Zinnecker 
&; Yorke 2007). In recent years, MIR shells or bubbles 



around the expanding Hll regions are recognised as the 
sites to observationally investigate the conditions of se- 
quential/triggered star formation ( Elmegreen &; Lada|1977 



Elmegreen 2010, and references therein) and the formation 



of new massive star(s) and/or cluster (s) as well. Two 
mechanisms have been proposed to explain the observed 
star formation due to influence of massive star(s): "collect 



and collapse" (Elmegreen & Lada 



1977 



Whitworth et al. 



1994) and radiation-driven implosion (RDI; Bertoldi 1989 
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Lefloch & Lazareff |1994| ). In the "collect and collapse" 
scenario, the Hll region expands and accumulates molecular 
material between the ionization and the shock fronts. With 
time the collected material becomes unstable and fragments 
into several clumps and lead to the formation of new 
generation of stars, as an effect of the shocks. In the RDI 
model, the expanding Hll region supply enough external 
pressure to initiate collapse of a pre-existing dense clump 
in the molecular material. 

In this work, we present a multi-wavelength study of 
a MIR bubble N14 associated with IRAS 1813 4-1652, from 

2007| around the 



the catalog of |Churchwell et al. ( 2006 



Galactic H 11 region G014.0-00. 1. The bubble N14 is situated 
at a near distance of 3.5 kpc ( |Beaumont &; Willia ms 2010) 
and is associated with the water maser near (~ 33 arcsec) 
to the IRAS position ( |Codella et al.|1994| ). [Lockman| ( [T989t 
reported the velocity of ionized gas (vlsr) to be about 
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36 km s 1 near to the IRAS position, using a hydrogen 



recombination line study. Beaumont &; Williams (2010) 



and Deharveng et al. (2010) studied several MIR bubbles 



including the N14 using James Clerk Maxwell Telescope 
(JCMT) 12 CO(J=3-2) line and APEX Telescope Large 
Area Survey of the Galactic plane at 870 /im (ATLASGAL) 
870 /im continuum observations, respectively. Beaumont & 
Williams ( |2010| also reported 20 cm Multi- Array Galactic 
Plane Imaging Survey (MAGPIS) radio continuum data 

~ fl2010j estimated 



around the N14. Beaumont & Williams 



molecular gas velocity to be about 40.3 km s~ associated 
with the bubble N14 with velocity dispersion of 2.8 km s _1 . 
They also listed about six 09.5 stars to produce observed 
MAGPIS 20 cm integrated flux (~ 2.41 Jy) for the Hll 
region associated with the bubble N14. Deharve ng et al.] 
(2010) found three dense clumps around the bubble N14 



and stated that the bubble is broken in the direction of the 
low density region. 

Previous studies on this region therefore clearly reveal 
the presence of molecular, cold dust as well as ionized 
emissions in the bubble. In this paper, we present multi- 
wavelength observations to study the interaction of Hn 
region with the surrounding interstellar medium (ISM). 
Our study will allow us to investigate the star formation 
especially the identification of embedded populations and 
also explore whether there is any evidence to form stars by 
the triggering effect of the Hll region around the bubble 
N14. 

In Section [2] we introduce the archival data and data 
reduction procedures used for the present study. In Sec- 
tion [3] we examine the structure of the MIR bubble N14 in 
different wavelengths and the interaction of massive stars 
with its environment using various Spitzer MIR ratio maps. 
In this section, we also describe the selection of young 
population, their distribution around the bubble, identi- 
fication of ionizing candidates and discuss the triggered 
star formation scenario on the borders of the bubble. In 
Section [4] we summarize our conclusions. 



2 AVAILABLE DATA AND DATA 
REDUCTION 

Archival deep near-infrared (NIR) HK S images and a cata- 
log around the bubble N14 were obtained from the UKIDSS 
6 th archival data release (UKIDSSDR6plus) of the Galac- 



tic Plane Survey (GPS) ( jLawrence et al.||2007| ). It is to be 
noted that there is no GPS J band observation available 
for the bubble N14. UKIDSS observations were made us- 



ing the UKIRT Wide Field Camera (WFCAM; ICasali et al. 



2007) and fluxes were calibrated using Two Micron All Sky 
Survey (2MASS; [Skrutskie et al.|2006) . The details of basic 
data reduction and calibration procedures are described in 
Dye et al. (2006) and Hodgkin et al. (2009), respectively. 



Magnitudes of bright stars (H ^ 11.5 mag and K s ^ 10.5 
mag) were obtained from the 2MASS, due to saturation of 
UKIDSS bright sources. Only those sources are selected for 
the study which have photometric magnitude error of 0.1 
and less in each band to ensure good photometric quality. 
We obtained narrow-band molecular hydrogen (H 2 ; 2.12 /im; 



1 - S(l)) imaging data from UWISH2 survey (Froebrich et 
al.|2011 ). We followed a procedure similar to that described 



by |Varrica tt (2011 ) to obtain the final continuum-subtracted 
H2 image using GPS K s image. 

The Spitzer Space Telescope Infrared Array Camera (IRAC 
(Chi (3.6 /im), Ch2 (4 .5 /mi), Ch3 (5.8 /im) and Ch4 (8.0 
/im) ; I Fazio et aL 2004 ) and Multiband Imaging Photometer 
(MIPS (24 /mi); Rieke et al.|2004| archival images were ob- 
tained around the N14 region from the "Galactic Legacy In- 
frared Mid-Plane Survey Extraordinaire" (GLIMPSE; |Ben- 



jamin et al.||2003| IChurchwell et al. 1120091) and "A 24 and 



70 Micron S urvey of the Inner Galactic Disk with MIPS" 
(MIPSGAL; |Carey et al.][2005| surveys. MIPSGAL 24 /im 
image is saturated close to the IRAS position inside the 
bubble N14. We used GLIMPSE-I Spring '07 highly reliable 
Point-Source Catalog and also performed aperture photom- 
etry on all the GLIMPSE images (plate scale of 0.6 arc- 
sec/pixel) using a 2.4 arcsec aperture and a sky annulus 
from 2.4 to 7.3 arcsec using; irafQ for those sources de- 
tected in images but the photometric magnitudes are not 
available in the GLIMPSE-I catalog. The IRAC/GLIMPSE 
photometry is calibrated using zero magnitudes including 
aperture corrections, 18.5931 (Chi), 18.0895 (Ch2), 17.4899 
(Ch3) and 16.6997 (Ch4), obtained from IRAC Instrument 
Handbook (Version 1.0, February 2010). 
We obtained 20 cm radio continuum map (resolution ~ 6 



arcsec) from Very Large Array MAGPIS survey (|Helfand et 



al. 2006| ) to trace the ionized region around the N14. The 
molecular 12 CO(J=3-2) (rest frequency 345.7959899 GHz) 
spectral line public processed archival data was also uti- 
lized in the present work. The CO observations (project 
id: M10BD02) were taken on 22 August 2010 at the 15 m 
JCMT using the HA RP array. Archival BOLOCAM 1.1 mm 
(Ag uirre et al.|20lT ) image (with effective FWHM Gaussian 
beam size of ~ 33 arcsec) was also used in the present work. 



3 RESULTS AND DISCUSSION 

3.1 Multi- wavelength view of the bubble N14 

Figure [l] shows the selected region (~ 12 x 8.6 arcmin 2 ) 
around the bubble N14, made of the 3-color composite 
image using GLIMPSE (8.0 /im (red) & 4.5 /im (green)) 
and UKIDSS K s (blue). The 8 /im band contains the two 
strongest PAH features at 7.7 /im and 8.6 /im, which are ex- 
cited in the photodissociation region (or photon-dominated 
region, or PDR). The PDRs are the interface between neu- 
tral & molecular hydrogen and traced by PAH emissions. 
PAH emission is also known to be the tracer of ioniza- 
tion fronts. Th e positions of IRAS 1 8134-1652 (+) and wa- 
ter maser (x) (Codella et al. 1994) are marked in the fig- 
ure. Figure [l] displays an almost circular morphology of the 
MIR bubble N14 prominently around the IRAS 18134-1652. 
These bubble structures are not seen in any of the UKIDSS 
NIR images, but prominently visible in all GLIMPSE im- 
ages. JCMT molecular gas (JCMT CO 3-2) emission con- 
tours are also overlaid on the Fig. [l] with 20, 40, 60, 80 and 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, USA 
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Figure 1. Three color composite image (size ~ 12 x 8.6 arcmin 2 ; central coordinates: 0:2000 = 18^ 16 m 18 s . 9, #2000 = -16 51 25 .2) 
of the selected region around the bubble N14, using Spitzer-GLIMPSF, images at 8.0 fim (red), 4.5 fim (green) and UKIDSS K s (blue) 
in log scale. Archival JCMT molecular 12 CO (3-2) gas emission contours are also overlaid on the image with 20, 40, 60, 80 and 95 % of 
the peak value i.e. 60.43 K km s -1 . The positions of the three 870 /im dust clumps from |Deharveng et al.| ( |2010| are marked by big red 
star symbols in the image. The scale bar on the top left shows a size of 1 pc at the distance of 3.5 kpc. The positions of IRAS 18134-1652 
(+) and water maser (x) are marked in the figure. 
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Figure 2. 20 cm contours in black color around the bubble N14 obtained from the MAGPIS survey, overlaid on a color composite image 
made using the 24 /im (red), 8 /im (green), and 3.6 /im (blue) images. The 20 cm contour levels are 40, 55, 70, 85 and 95 % of the peak 
value i.e. 0.0115 Jy/ beam. BOLOCAM 1.1 mm emission is also shown by solid yellow contours with 20, 30, 40, 55, 70, 85 and 95 % 
of the peak value i.e. 1.07 Jy/ beam. The peak positions of BOLOCAM 1.1 mm clumps are shown by triangle symbols obtained from 
the Bolocam Galactic Plane Survey (BGPS) catalog. The other marked symbols are similar to those shown in Fig. ^ (see text for more 
details). MIPS 24 /im image is saturated near the IRAS position. 
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95 % of the peak value i.e. 60.43 K km s 1 . The peak posi- 



line between 39 to 41.5 km s . These values are also consis- 



tions of the three detected 870 /im dense clumps (Deharveng 
et al.|20To| are also marked by a big star symbol in the fig- 
ure. Figure [2] shows a color composite image made using 
MIPSGAL 24 /im (red), GLIMPSE 8 /im (green), and 3.6 
/im (blue) images, overlaid by BOLOCAM 1.1 mm emission 
by solid yellow contours. MIPS 24 /im image is saturated 
near to the IRAS position, but a few point sources are also 
seen around the bubble. MAGPIS 20 cm radio continuum 
emission is also overlaid in Figure [2] by black contours. The 
cold dust emission at 1 . 1 mm is very dense and prominent at 
the borders of the bubble. It is clearly seen that the peaks 
of 870 /im, 1.1 mm dust emission and CO molecular gas 
emissions are spatially coincident along the borders of the 
bubble (see Figs.[l]&|2}. The 24 /xm and 20 cm images trace 
the warm dust and ionized gas in the region, respectively. It 
is obvious from Figure [2] that the PDR region (traced by 8 
/xm) encloses 24 /im dust and 20 cm ionized emissions inside 
the bubble and indicates the presence of dust and gas in and 
around the Hll region (see e. g. Wat son et al.|2008 for N10, 
N21, and N49 bubbles). 

3.2 PAH emission and the Collected material 

In recent years, Spitzer-IRAC bands and ratio maps are 
utilized to study the interaction of massive stars with its 
immediate environment (Povich et al 2007). The IRAC 



bands contain a number of prominent atomic and molec- 
ular lines/features such as H2 lines in all channels (see Ta- 
ble 1 from |Smith fe Rosen|[2U05t , Bra 4.05 /im (Ch2), Fell 



5.34 /im (Ch3), A m 6.99 /im and Ann 8.99 /im (Ch4) (see 
Reach et al. 2006). We know that the Spitzer-IRAC bands, 
Chi, Ch3 and Ch4, contain the PAH features at 3.3, 6.2, 
7.7 and 8.6 /im, whereas Ch2 (4.5 /im) does not include any 
PAH features. Therefore, IRAC ratio (Ch4/Ch2, Ch3/Ch2 
and Chl/Ch2) maps are being used to trace out the PAH 



features in massive star forming (MSF) regions (e.g. [Povich 


et al 2007; Watson et 


al.||2008| |Kumar & Anandarao 2010; 


Dewangan et al. 


2012 


) due to UV radiation from massive 



star(s). In order to make the ratio maps, we generate resid- 
ual frames for each band removing point sources by choosing 
an extended aperture (12.2 arcsec) and a larger sky annulus 
(14.6 - 24 .4 arcsec; |Reach et aLpOO^ in IRAF/DAOPHOT 



software (Stetson 1987). These residual frames are then sub- 
jected to median filtering with a width of 4 pixels and 
smoothing by 9 x 9 pixels using the boxcar algorithm. Fig- 
ures [3^1 and [3j) represent the IRAC ratio maps, Ch3/Ch2 
and Ch4/Ch2 around the bubble N14, respectively. The ra- 
tio contours are also overlaid on the maps for better clarity 
and insight (see Fig. [3]). Both the ratio maps clearly trace 
the prominent PAH emissions and subsequently the extent 
of PDRs in the region. IRAC ratio maps reveal that the 
bubble is traced by the PAH emission following an almost 
circular morphology except in south-west direction towards 
the low molecular density environment, which was also re- 



ported by Deharveng et al. (2010) 



The emission contours of dust (ATLASGAL 870 /im and 
BOLOCAM 1.1 mm) and molecular gas (JCMT CO 3-2) 
exhibit the evidence of collected material along the bubble 
(see Figs, [I] [2] and [3} . |Deharveng et al.| fl2010| ) tabulated 
the positions of three 870 /im dust condensations with their 
respective molecular velocity using the NH 3 (1,1) inversion 



tent with Beaumont & Williams (2010), who also reported 
velocity for the bubble N14 (~ 40.3 km s _1 ). These velocity 
ranges of molecular gas are also compatible with the ion- 
ized gas velocity (~ 36 km s _1 ) obtained by the hydrogen 



recombination line study around the Hn region (Lockman 



1989), which confirms the physical association of the molec- 
ular material and the H II region. It is also to be noted that 
the peaks of cold dust and molecular gas emission at differ- 
ent locations possibly indicate the fragmentation of collected 
materials into different individual clumps around the bub- 
ble. The evidence of collected material along the bubble is 
further confirmed by the detection of the H2 emission (see 
Figure [4). Figure [4] represents the continuum-subtracted H2 
image at 2.12 /im and reveals that the H2 emission sur- 
rounds the Hn region along the bubble, forming a PDR 
region, which may be collected due to the shock. In brief, 
the PAH emission, cold dust emission, molecular CO gas 
and shocked H2 emissions are coincident along the bubble. 



3.3 Photometric analysis of point-like sources 
towards N14 

In order to trace ongoing star formation activity around 
the bubble N14, we have identified YSOs using NIR and 
GLIMPSE data. 



3.3.1 Selection of YSOs 
We have used 



Gutermuth et al. (2009) criteria based 



on four IRAC bands to identify YSOs and various pos- 
sible contaminants (e.g. broad-line AGNs, PAH-emitting 
galaxies, shocked emission blobs/knots and PAH-emission- 
contaminated apertures). These YSOs are further classified 
into different evolutionary stages (i.e. Class I, Class II, Class 
III and photospheres) using slopes of the IRAC spectral 
energy distribution (SED). Fig. |5^l shows the IRAC color- 
color ([3.6]-[4.5] vs [5.8]-[8.0]) diagram for all the identified 
sources. We find 33 YSOs (15 Class 0/1; 18 Class II), 621 
photospheres and 78 contaminants in the selected region 
around the bubble N14. The details of YSO classifications 
can be found in Dewangan & Anandarao ( 2011| . We have 
also applied criteria ([3 . 6]-[4.5] 0.7 and [4.5 ]-[5.8] = 0.7; 
jHartmann et al.| ( |2005| ); jGetman et al.| ( j2007| )) to identify 
protostars (Class I) among the sources, which are detected 
in three IRAC/GLIMPSE bands, but not in the 8.0 /im band 
and rest of the remaining sources are subjected to SED crite- 



ria (see |Dewangan &; Anandarao|20lT ) to select the Class II 
and Class III sources. We identify 186 additional YSOs (10 
Class 0/1; 176 Class II) through color-color diagram using 
three GLIMPSE bands in the region (see Fig. [5)3). 
GLIMPSE 3.6 and 4.5 /im bands are more sensitive for point 
sources than GLIMPSE 5.8 and 8.0 /im images. Therefore, 
a larger number of YSOs can be identified using a combina- 
tion of UKIDSS NIR HK S photometry with GLIMPSE 3.6 
and 4.5 /im (i.e. NIR-IRAC) photometry, where sources are 
not detected in IRAC 5.8 and/or 8.0 /im band (Gutermuth 
et al.l 2009 ) . We followed the criteria given by 



Gutermuth 



et al.| ( |2009| ) to identify YSOs using H, K s , 3.6 and 4.5 /im 
data. We have found 199 additional YSOs (185 Class II and 
14 Class I) using NIR-IRAC data (see Fig. ^jp). 
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Figure 3. a) GLIMPSE Ch3/Ch2 ratio map of the region around N14 (similar area as shown in Fig. [TJ. The ratio Ch3/Ch2 value is 
found to be "7-8" and "8.5-10" in the bubble interior and for the brightest part of the PDR respectively. The Ch3/Ch2 ratio contours 
are also overlaid on the image with a level of 8.3, a representative value between "7-8.5". b) The Ch4/Ch2 ratio map of the region 
is shown here. The ratio Ch4/Ch2 value for the brightest part of the PDR and the interior is "30-34" and "23-25" respectively. The 
Ch4/Ch2 ratio contours are also overlaid on the image with a level of 27.8, a representative value between "25-30". The white dashed 
box is shown as a zoomed-in view in Fig. ^] MAGPIS 20 cm radio contours in yellow color are also overlaid in both the figures. The 
marked symbols and levels of radio contours on the figure are similar to those shown in Fig. [T] 



Finally, we have obtained a total of 418 YSOs (379 Class II 
and 39 Class I) using NIR and GLIMPSE data in the region. 



3.3.2 Spatial distribution of YSOs 

To study the spatial distribution of YSOs, we generated the 
surface density map of all YSOs using a 5 arcsec grid size, 
following the same procedure as given in [Gutermuth et all] 
(2009). The surface density map of YSOs is constructed 



using 6 nearest-neighbor (NN) YSOs for each grid point. 



Fig. [5]i shows the spatial distribution of all identified YSOs 
(Class I and Class II) in the region. The contours of YSO 
surface density are also overlaid on the map. The levels 
of YSO surface density contours are 6 (2.1cr), 9 (3.2a), 14 
(4.9a), 20 (7.0a) and 26 (9.1a) YSOs/pc 2 , increasing from 
the outer to the inner region. We have also calculated the 
empirical cumulative distribution (ECD) as a function of 
NN distance to identify the clustered YSOs in the region. 
Using the ECD, we estimate the distance of inflection d c 
= 0.54 pc (0.0088 degrees at 3.5 kpc) for the region for a 
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Figure 4. Zoomed-in inverted gray-scale continuum-subtracted H2 image at 2.12 /im around N14 (size ~ 7.2 x 5.7 arcmin 2 ; as shown 
by a dashed box in Fig.[3j)). The H2 image clearly indicates the detected H2 emission (dark features as seen in GLIMPSE images) in the 
selected region around N14. The marked symbols are similar to those shown in Fig. ^ The continuum-subtracted H2 image is processed 
through median filtering with a width of 4 pixels and smoothened by 4 x 4 pixels using boxcar algorithm to trace the faint features in 
the image. 



surface density of 6 YSOs/pc 2 (see Dewangan &; Anandarao 
2011 for details of d c and ECD). We find that about 23% 



our selected sources from Robitaille et al. (2008), wherever 



(98 out of 418) YSOs are present in clusters. Figure [5]i 
reveals that the distribution of YSOs is mostly concentrated 
in and around the bubble, having peak density of about 
20 YSOs/pc 2 (see also Figure [6), while YSO density of 
about 8 YSOs/pc 2 (2.8<j) is also seen around the east and 
west regions close to the low density molecular gas and 
dust emission (see Figure [6]). It is also to be noted that 
the YSOs clustering is seen along the PDR region, on 
the borders of the bubble. The correlation of cold dust, 
molecular gas, ionized gas and YSO surface density is shown 
in Figure [6] The association of YSOs with the collected 
materials around the region further reveals the ongoing star 
formation on the borders of the bubble. 

In addition, we have checked the possibility of intrin- 
sically "red sources" contamination, such as asymptotic 
giant branch (AGB) stars in our YSO sample. Recently, 
Robitaille et al.| ( |2008| prepared an extensive catalog of such 
red sources based on the Spitzer GLIMPSE and MIPSGAL 
surveys. They showed that two classes of sources are well 
separated in the [8.0 - 24.0] color space such that YSOs are 
redder than AGB stars in this space (see also |Whitney et 
ah][2008). Firstly, we applied red source criteria for all our 
selected YSOs having 4.5 /im and 8.0 /im detections. We 
find 19 out of 33 YSOs (~ 57%) as possible intrinsically 
"red sources" and then further utilized [8.0 - 24.0] color 
space to identify the AGB contaminations from selected 
"red sources". We have taken MIPS 24 /im magnitude for 



it is available and also extracted a few more sources from 
archival MIPSGAL 24 /im image. We performed aperture 
photometry on MIPSGAL 24 /im image using IRAF with 
a 7 arcsec aperture & a sky annulus from 7 to 13 arcsec. 
Zero points and aperture corrections were adopted from 
MIPS Instrument Handbook Version 3, March 2011, for 
selected aperture. MIPS 24 /im image is saturated close to 
the IRAS position, however, a few point sources are seen 
around the bubble. Therefore, we obtained MIPS 24 /im 
magnitudes for only 9 YSOs (7 Class I and 2 Class II) iden- 
tified as red sources. Finally, we identified two likely AGB 
contaminations out of 9 red sources, following the criteria 
suggested by ||Robitaille et alT (2008), which are located 
away from the identified YSO clusters (see Figure |6|. We 
have not considered these likely AGB contaminations for 
further analysis. 



3.3.3 SED modeling of Class I YSOs 

In this subsection, we present SED modeling of all 37 

identified Class I YSOs (designated as si, ,s37) as well 

as two selected YSOs (s38 and s39) associated with the 
peak of molecular gas and dust clumps around the bubble, 
in our selected region around N14, to derive their various 
physical parameters using an on-line SED modeling tool 
(Robitaille et al. 112006 20071). It is interesting to note that 



"s39" is a deeply embedded source (prominently seen in 24 
/im image), detected only in 5.8 /im and longer wavelength 
bands. NIR and Spitzer IRAC/GLIMPSE photometric 
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Figure 5. a) Color-color diagram (CC-D) using Spitzer-IUAC four bands for all the sources identified within the region shown in Fig.^ 
The extinction vector for = 5 mag is shown by an arrow, using average extinction law from |Flaherty et al.| j2007| ). The dots in gray 
color around the centre (0,0) locate the stars with only photospheric emissions. The open squares (black), open triangles (blue) and 
open circles (red) represent Class III, Class II and Class 0/1 sources respectively, classified using the olirac criteria. The "x" symbols 
in magenta color show the identified PAH-emission-contaminated apertures in the region. The PAH galaxy contamination source is also 
marked by a violet cross (+) symbol in the diagram (see the text), b) CC-D of the sources detected in three IRAC/ GLIMPSE bands, 
except 8.0 /im. The small filled cirlces (black), open triangles (blue) and open circles (red) represent Class III, Class II and Class 0/1 
sources respectively (see the text for YSOs selection critria). c) Figure shows the de-reddened [K s - [3.6]] o vs [[3.6] - [4.5]] o CC-D using 
NIR and GLIMPSE data. The selected region shown by the solid lines represents the location of YSOs. The extinction vector for 
= 2 mag is shown by an arrow, calculated using the average extinction law from Flah erty et al.| ( |2007| >. Open red circles and open 
blue triangles represent Class I and Class II sources respectively, d) Plot represents the spatial distribution of all identified YSOs in the 
selected region around N14. The YSO surface density contours are plotted for 6, 9, 14, 20 and 26 YSOs/pc 2 , from outer to inner side 
(see text for details). The open circles and open triangles show the Class I and Class II sources respectively. The YSOs identified using 
four IRAC, three IRAC and NIR-IRAC data are shown by red, black and blue colors respectively. The JCMT 12 CO (3-2) contours in 
gray color and marked symbols are similar to those shown in Fig. ^ 
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Figure 6. The spatial distribution of YSOs, dust emission and molecular gas in the selected region around N14. BOLOCAM 1.1 mm 
image is overlaid with contours of YSO surface density (blue; similar levels as presented in Fig. |5jl) and JCMT CO 3-2 (yellow; similar 
levels as shown in Fig. [1} data. All Class I YSOs (see subsection |3.3.3| are also marked (red circles) and labeled as 1,....,37 on the image 
(see Table [TJ. Two selected sources (s38 and s39) are shown by filled yellow squares and labeled, which are associated with the peaks of 
molecular gas and dust emission around the bubble. The position of two likely AGB contaminations is also marked by orange asterisks 
on the image. Two boxes are shown as a zoomed-in view in Fig. [7] 



magnitudes for these selected YSOs are listed in Table [T] 
along with IRAC spectral indices (q>irac) and are also 
labeled in Fig. [6] Spitzer 24 /xm magnitudes are also listed 
for selected YSOs, wherever it is available. Figs. [T^l and 
[TJ) show the zoomed-in 3 color composite image using 
GLIMPSE (5.8 fim (red) & 4.5 /mi (green)) and UKIDSS 
K s (blue) around the east of the bubble close to the peak of 
a dense clump and around the bubble, respectively. Fig. [T^l 
clearly exhibits the location of a deeply embedded source 
"s39" along with other identified Class I sources close to 
the peak of a cold dust emission. Fig. [7}) presents the 
zoomed-in view around the bubble with the positions of a 
few identified YSOs (such as sl7, sl8 and s38, along with 
other identified Class I YSOs). 

IRAC spectral indices were calculated using a least squares 
fit to the IRAC flux points in a log(A) versus log(AFA) 
diagram for those sources that are detected in atleast 
3 IRAC bands (see Dewangan & Anandarao 2011[ for 



details). The SED model tool requires a minimum of three 
data points with good quality as well as the distance to the 
source and visual extinction value. These models assume 
an accretion scenario with a central source associated 
with rotationally flattened infalling envelope, bipolar 
cavities, and a flared accretion disk, all under radiative 
equilibrium. The model grid consists of 20,000 models 
of two-dimensional Monte Carlo simulations of radiation 
transfer with 10 inclination angles, resulting in a total of 
200,000 SED models. The grid of SED models covers the 
mass range from 0.1 to 50 M©. Only those models are 



shown in Fig. [8] for 9 out of 39 selected YSOs associated 
with the molecular and dust clumps. The weighted mean 
values of the physical parameters (mass and luminosity) 
along with the standard deviations derived from the SED 
modeling for all the selected sources are given in Table [I] 
The SED results clearly show that the YSOs having higher 
luminosity represent more massive candidates. The table 
also contains the number of models that satisfy the x 2 
criterion as mentioned above. The derived SED model 
parameters show that the average values of mass and 
luminosity of the 39 selected sources are about 4.8 M© and 
1817.1 L©, respectively. Our SED result shows that the 
source "s4" is the most luminous and massive YSO (~ 20.5 
M©) among all selected YSOs away from the bubble and 
is saturated in the GLIMPSE 8 /xm image. It is however 
tabulated as an OH selected AGB/post-AGB candidate 



selected that satisfy the criterion x 



"Kbest 



< 3, where % 



is taken per data point. The plots of SED fitted models are 



by |Sevenster| (2002) using 1612 MHz masing OH line 
profile. The YSO surface density contours clearly trace a 
clustering of Class I YSOs (s27-28, s30-33 and s39) with ~ 
20 YSOs/pc 2 associated with the dense dust clump at the 
eastern border of the bubble N14 along with other peaks of 
YSO surface density (see Figs|6][7^ and Table [lj. It is found 
that about 5 young massive embedded YSOs (s6, sl8, s23, 
s28 and s31) with a mass range of 8 to 10 M© and about 
15 intermediate mass YSOs (s7, s9-10, sl2, sl4, sl6-17, 
sl9, s22, s25-26, s30, s33, s38-39) with mass range of 3 to 7 
M© are associated with the molecular and dust fragmented 
clumps at the borders of the bubble (see Figs. [6] [7] and 
Table [TJ. It is interesting to note that the sources sl8, 
s28, s31, s38 and s39 are associated with the peak of dust 
clumps at the border of the bubble and three of them 
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Right Ascension (J2000) 

Figure 7. Three color composite image (UKIRT K s (blue) and GLIMPSE 4.5 jum (green), 5.8 /im (red)) of the two zoomed-in regions 
around the bubble N14. a) Figure clearly exhibits the location of a deeply embedded source "s39" along with other identified Class I 
sources close to the peak of cold dust emission (see white solid box in Fi g. [6| . b) The MAGPIS 20 cm contours are overlaid on the 
zoomed-in region (dashed box in Fig. [6} with similar levels as shown in Fig. |2[close to the bubble N14 (size ~ 3.7 x 3.1 arcmin 2 ). The 8 
selected main sequence stars as probable ionizing candidates (see subsection |3. 3. 4| ) are also marked (black star symbols) and labeled as 
1,....,8 on the image (see Table [2j. The diamonds in red color are identified as the O-type ionizing stars (#4 and #6) inside the bubble. 
Some of Class I YSOs in yellow color are also marked and labeled on the figure. 



(sl8, s28 and s31) are possibly young massive protostars. 
Finally, the SED modeling results favor ongoing star 
formation around the region with detection of YSOs as well 
as some massive protostars in their early phase of formation. 



Recently, Kryukova et al. (2012) studied a relation- 



ship between bolometric luminosity and MIR luminosity 



(integrated from 1.25 /im to 24 /im) of Spitzer identified 
protostars in nine nearby molecular clouds, independent of 
SED modeling. We have also computed the MIR luminosity 
(Lmih) of our selected sources from integrating the SED 
using their available photometric infrared data (see Table [I] 
for Lmir)- However, the estimated MIR luminosity of our 
selected sources is underestimated because of the lack of 
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Figure 8. The SED plots of 9 out of 39 selected YSOs associated with the molecular and dust clumps are shown here. Filled circles are 
observed fluxes of good quality (with filled triangles as upper limits) taken from the archives or published literature (see text for more 
details) and the curves show the fitted model with criteria x 2 - Xbest < ^. The thin black curve corresponds to the best fitting model. 
The dashed curves represent photospheric contributions. 



24 fim and longer wavelength data for most of the sources. 
Therefore, we prefer to use model derived SED physical pa- 
rameters (like mass and luminosity) over luminosity derived 
independently of SED modeling for our selected sources. It 
is known that the physical parameters obtained from SED 
modeling are not unique but indicative. We therefore used 
only mass of selected sources for their relative comparison 
in terms of mass. 



3.3.4 Ionizing Candidates 

The presence of both PAH emission around the bubble and 
radio emission inside the bubble clearly suggest the pres- 
ence of the ionizing source(s) with UV radiation close to 
the centre of the bubble. One can find more details regard- 
ing identification of ionizing candidate (s) inside the bubble 
Pomares et~aT1 ( |2QQ9| ) and |Ji et al.| fl2012| ). It is to be 



we followed suggestions of Pomares et al. (2009) and Ji et 



noted that no ionizing stars are reported for this bubble. 



Beaumont &; Williams (2010) listed about six 09.5 stars 
to produce observed 20 cm integrated flux (~ 2.41 Jy) for 
the Hn region associated with the bubble N14. Therefore, 



al. (2012) to trace the ionizing candidate (s) inside the bub- 



ble N14 using 2MASS and GLIMPSE photometric magni- 
tudes. We selected 8 candidates (designated as #1, ,#8) 

to search for O-type star(s) inside the bubble based on their 
detections in J band to 5.8 /xm or longer wavelength bands 
(see Fig.[7]3 and Tablefg]). All selected candidates are marked 
and labeled in Fig. ^p. We calculated the absolute JHK S 
magnitude for each candidate using their 2MASS apparent 
J, H, and K s magnitudes. 2MASS photometry is used here 
due to non availability of GPS J band photometry. We used 
a distance of 3.5 kpc and estimated the extinction for each 
source from the NIR color color diagram (CC-D). We fol- 
lowed the extinction law given by Indeb etouw et al.] ([2005 ) 
(Aj/Av = 0.284, Ah/Av = 0.176 and A K /A V = 0.114) 
and used the intrinsic colors (J - H)o (= -0.11) and (H - 
K) (= -0.10) obtained from |Martins fe Plezj p006| . We 
compared the derived absolute JHK S magnitudes with those 
listed by Martins & Plez (2006) for our selected candidates, 



and found two O-type candidates (#4 and #6) inside the 
bubble. We also checked their evolutionary stages using CC- 
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D (see subsection 3.3.1) and found that all the sources are 
main sequence stars except source #3, which is identified as 
a Class I YSO and designated as sl4 in Table [I] The posi- 
tions of the 8 selected candidates are tabulated in Table [2] 
with their 2MASS NIR & GLIMPSE apparent magnitudes, 
calculated visual extinction, estimated absolute JHK S mag- 
nitudes and the possible spectral class obtained from the 
comparison of absolute magnitude of each candidate with 



the listed values of |Martins & Plez| |2QQ6) 



3.4 Star formation scenario 

We have found evidence of collected material along the bub- 
ble and also ongoing formation of YSOs on the borders of 
the bubble. The YSO clusters and embedded YSOs discov- 
ered are associated with the peak of molecular and cold dust 
material collected on the borders of the bubble. The mor- 
phology and distribution of YSOs suggest that the bubble 
N14 is a site of star formation possibly triggered by the ex- 
pansion of the H II region. In recent years, the triggered star 
formation process, especially "collect and collapse" mecha- 
nism has been studied extensively on the borders of many 
Hn regions such as Sh 2-104, RCW 79, Sh 2-212, RCW 120, 



Sh 2-217, G8.14+0.23 (|Deharveng et al.||2003 


2008, |2009| 


Zavagno et al.||2006| |2010| |Brand et al.||2011||] 


Dewangan et 


al.|2012|. In order to check the "collect and collapse" process 



as the triggering mechanism around N14, we have calculated 
the dynamical age (tdyn) of the Hll region and compared it 



with an analytical model by Whitworth et al. (1994). We 



have estimated the age of the Hll region at a given radius 
R, using the following equation (Dyson & Williams 1980| : 



tdyn 



4Rs 
7c 3 



7/4 



(1) 



where c s is the isothermal sound velocity in the ionized gas 
(c s = 10 km s" 1 ) and R s is the radius of the Stromgren 
sphere, given by R s = (3 /47m 2 ) a^) 1 7/3 , where the ra- 
diative recombination coefficient a B = 2.6 x 10" 13 (10 4 
K/T) 7 cm 3 s" 1 (Kwan|l997k. In this calculation, we have 
used as — 2.6 x 10" 13 cm^ s" 1 for the temperature of 10 4 
K. N uv is the total number of ionizing photons per unit 
time emitted by ionizing stars and "no" is the initial par- 
ticle number density of the ambient neutral gas. We have 
adopted the Lyman continuum photon flux v alue (N uv — 
of 2.34 x 10 48 ph s" 1 (logN u 



48.36) from Beaumont & 



Williams (2010) for an electron temperature, distance and 
integrated 20 cm (1.499 GHz) flux density of 10 4 K, 3.5 kpc, 
and 2.41 Jy respectively. We calculated the mean H2 num- 
ber density near the Hll region using archival 12 CO zeroth 
moment map or column density map. In general, the zeroth 
moment map (in unit of K km s" 1 ) is created by integrat- 
ing the brightness temperature over some velocity range. We 
derived the column density using the formula A^h 2 (cm -2 ) 
= X x Wco, assuming the molecular clumps/cores are ap- 
proximately spherical in shape. We used the CO-H 2 conver- 
sion factor (also called X factor) for dense gas as 6 x 10 20 



cm" 2 K km" 1 s from |Shetty et al.|fl201l| ). We find almost 
a closed circular structure of integrated CO line emission 
(Wco) (typical value of ~ 23.43 K km s" 1 ) using archival 
column density map, which traces well the bubble around 
the Hn region. The mean H2 number density is obtained 



to be 2071.5 cm" 3 using the relation Nn 2 (cm" 2 )/L (cm), 
where L is the molecular core size of about 6.79 x 10 18 cm (~ 
2.2 pc) near the Hll region. The estimated mean H2 number 
density could be under-estimated because of not considering 
the fact that the 12 CO (3-2) transition can be (partly) op- 
tical thick in nature and also assuming a spherical structure 
in the calculation. Using N uu , a radius of the Hll region (R 
=) 2.48 pc, and n = 2071.5 cm" 3 , we have obtained tdyn 
^0.74 Myr using Equation 1. The timescale of Hll region 
expansion is actually not an easy issue. There is the long- 
standing problem of the lifetime of Galactic ultra-compact 
(UC) H 11 regions. |Wood & ChurchweH] ( |1989| ) observation- 
ally reported that the typical lifetime of UC H II regions is 
about 10 5 years, but the lifetime of an expanding UC Hll 
region (i.e. the timescale to reach pressure equilibrium with 
the surrounding environment) is estimated to be about 10 4 
years, for 10 km s" 1 expansion velocity of ionized gas, with 
a typical radius of ~ 0.1 pc. So, the difference between these 
timescales is known as the lifetime problem of UC H 11 re- 
gions. This lifetime problem might have implications in the 
present case even if it is only the UC H II region which is liv- 
ing too long. Therefore, our estimated timescale should be 
considered with a caution. Following, the Whitw orth et al.| 
( 1994 ) analytical model for the "collect and collapse" pro- 
cess, we have estimated a fragmentation time scale (tf rag ) of 
1.27 - 2.56 Myr for a turbulent velocity (a s ) of 0.2 - 0.6 km 
s" 1 in the collected layer. We find that the dynamical age is 
smaller than the fragmentation time scale for no = 2071.5 
cm" 3 . We have checked the variation of tfrag and td yn with 
initial density (no) of the ambient neutral medium and found 
that if tdyn is larger than tf rag , then ambient density (no) 
should be larger than 3610, 5710, 6700 and 7700 cm" 3 for 
different turbulent velocity (a s ) values of 0.2, 0.4, 0.5 and 0.6 
km s" 1 respectively. We have also estimated the kinemat- 
ical time scale of the molecular bubble of about 1.57 Myr 
(~ 4.5 pc/ 2.8 km s" 1 ), assuming bubble size of about 4.5 
pc and velocity dispersion ~ 2.8 km s" 1 from 12 CO(J=3-2) 
map (see |Beaumont fe W illiams 2010). The comparison of 
the dynamical age of the Hll region with the kinematical 
time scale of expanding bubble and the fragmentation time 
scale, does not support the fragmentation of the molecu- 
lar materials into clumps due to "collect and collapse" pro- 
cess around the bubble. Also, the average age of Class 0/1 



sources is reported to be about 0.10-0.44 Myr (see | Evans et| 
|al.|2009| , which is less than the fragmentation time scale of 
the molecular materials into clumps. Therefore, we suggest 
the possibility of triggered star formation by compression 
of the pre-existing dense clumps by the shock wave and/or 
small scale Jeans gravitational instabilities in the collected 
materials. The YSO surface density contours clearly trace 
a clustering of Class I YSOs (s27-28, s30-33 and s39), with 
three intermediate mass (s30, s33 and s39) and two massive 
embedded sources (s28 and s31) associated with the dense 
dust clump at the eastern border of the bubble N14 (see 
Fig.[7ji and Table [TJ. We have also found that the sl8 & s38 
sources are also associated with the peak of different dust 
clumps around the bubble and the source "sl8" is identified 
as a new young massive protostar. 
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4 CONCLUSIONS 

We have explored the triggered star formation scenario 
around the bubble N14 and its associated Hn region us- 
ing mult i- wavelength observations. We find that there is 
a clear evidence of collected material (molecular and cold 
dust) along the bubble around the N14 region. The surface 
density of YSOs reveals ongoing star formation and clus- 
tering of YSOs associated with the borders of the bubble. 
We conclude that the YSOs are being formed on the border 
of the bubble possibly by the expansion of the Hn region. 
We further investigated the "collect-and-collapse" process 
for triggered star formation around N14 using the analyti- 



cal model of Whitworth et al. (1994). We have found that 



the dynamical age (~ 0.74 Myr) of the H n region is smaller, 
and the kinematical time scale of the bubble (~ 1.57 Myr) 
is comparable to the fragmentation time scale (~ 1.27 - 2.56 
Myr) of accumulated gas layers in the region for 2071.5 cm -3 
ambient density. The comparison of the dynamical age with 
the kinematical time scale of expanding bubble and the frag- 
mentation time scale does not support the fragmentation 
of the molecular materials into clumps due to the "collect 
and collapse" process around N14, but suggests the possi- 
bility of triggered star formation by compression of the pre- 
existing dense clumps by the shock wave and/or small scale 
Jeans gravitational instabilities in the collected materials. 
The YSO surface density contours clearly trace a clustering 
of Class I YSOs (~ 7 Class I sources with - 20 YSOs/pc 2 ) 
associated with the dense dust clump at the eastern border 
of the bubble N14. Also 5 young massive embedded proto- 
stars (about 8 to 10 Mq) and 15 intermediate mass (about 
3 to 7 M©) Class I YSOs are associated with the dust and 
molecular fragmented clumps at the borders of the bubble. 
It seems that the expansion of the H 11 region is also leading 
to the formation of these intermediate and massive Class I 
YSOs around the bubble N14. 
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Table 1. NIR, Spitzer IRAC/GLIMPSE and MIPS 24 /xm photometric magnitudes are listed for all selected YSOs in the region 
around the bubble (see the text). Spectral index (mirac) of sources was obtained having atleast 3 IRAC wavelengths magnitude. MIR 
luminosity (Lmir) of selected sources was obtained by integrating the SED using the available photometric infrared data. Weighted 
mean of luminosity (L*) and mass (M*) along with no. of models of all selected sources are also listed in the table, which are obtained 
by SED modeling (see the text). 
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Table 2. NIR and Spitzer IRAC/GLIMPSE photometric magnitudes are listed for 8 selected stars identified as probable ionizing 
candidate(s) inside the bubble N14. The estimated visual extinction Ay and absolute JHK S magnitudes for each source is also tabulated. 



ID 


RA 


Dec 


J 


H 


K s 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


Ay 


Mj 


M H 




O-type star 




[2000] 


[2000] 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 










1 


18:16:18.79 


-16:50:55.2 


14.78 


13.30 


12.68 


12.01 


11.72 


10.12 




13.14 


-1.67 


-1.73 


-1.53 




2 


18:16:22.14 


-16:50:19.8 


13.24 


12.02 


11.45 


11.03 


10.89 


10.45 




11.52 


-2.75 


-2.73 


-2.58 
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18:16:22.50 


-16:51:05.5 


13.07 


11.84 


11.32 


10.77 


10.43 


9.68 


7.39 


11.15 


-2.82 


-2.84 


-2.67 
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18:16:19.98 


-16:50:57.3 


12.64 


10.90 


10.09 


9.47 


9.32 


9.39 




15.89 


-4.59 


-4.62 


-4.44 


05V-04V 
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18:16:16.86 


-16:50:34.2 


15.95 
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9.45 


7.95 


7.84 


7.40 


7.61 


38.69 


-7.76 


-7.94 


-7.67 
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12.05 
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10.33 


9.84 


9.71 


9.62 


8.21 


11.19 


-3.85 


-3.79 


-3.66 


08V-07.5V 
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9.31 


7.77 


7.08 


6.83 


6.71 


6.51 


6.54 


13.98 
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-7.41 


-7.23 
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18:16:21.50 


-16:51:44.9 


15.35 


14.33 


13.80 


11.68 


11.20 


9.57 




10.27 


-0.29 


-0.20 


-0.09 





2012, ApJ, 144, 31 
Kwan J., 1997, ApJ, 489, 284 

Lawrence A., Warren S. J., Almaini O., et al. 2007, MN- 

RAS, 379, 1599 
Lefloch, B., Lazareff, B. 1994, A&A, 289, 559 
Lockman F. J., 1989, ApJS, 71, 469 
Martins F., Plez B., 2006, A&A, 457, 637 
Pomares M., Zavagno A., Deharveng L., et al., 2009, A&A, 

494, 987 

Povich M. S., Stone J. M., Churchwell E., et al., 2007, ApJ, 
660, 346 

Reach W. T., Megeath S. T., Cohen M., et al. 2005, PASP, 
117, 978 

Reach W. T., Rho J., Tappe A., et al., 2006, AJ, 131, 1479 

Rieke G. H. et al., 2004, ApJS, 154, 25 

Robitaille T. P., Whitney B. A., Indebetouw R., Wood K., 



Denzmore P., 2006, ApJS, 167, 256 
Robitaille T. P., Whitney B. A., Indebetouw R., Wood K., 

2007, ApJS, 169, 328 
Robitaille T. P., Meade M. R., Babler B. L., et al., 2008, 

AJ, 136, 2413 
Sevenster M. N., 2002, ApJ, 123, 2772 
Shetty, R., Glover, S. C, Dullemond, C. P., Klessen, R. S., 

2011, MNRAS, 412, 1686 
Skrutskie M. F., Cutri R. M., Stiening R., Weinberg M. D., 

Schneider S., Carpenter J. M., Beichman C, et al., 2006, 

AJ, 131, 1163 
Smith M. D., Rosen A., 2005, MNRAS, 357, 1370 
Stetson P. B., 1987, PASP, 99, 191 
Varricatt W. P., 2011, A&A, 527, 97 

Watson C, Povich M. S., Churchwell E. B., et al., 2008, 
ApJ, 681, 1341 



14 L. K. Dewangan & D. K. Ojha 



Whitworth A. P., Bhattal A. S., Chapman S. J., Disney 
M. J., Turner J. A., 1994b, MNRAS, 268, 291 

Whitney B. A., Sewilo M., Indebetouw R., et ah, 2008, AJ, 
136, 18 

Wood, D. O. S., Churchwell, E., 1989, ApJS, 69, 831 
Zavagno A., Deharveng L., Comeron F., et al., 2006, A&A, 
446, 171 

Zavagno A., Russeil D., Motte F., et al., 2010, A&A, 518, 
L81 

Zinnecker H., Yorke H. W., 2007, ARA&A, 45, 481 



